arXiv:1501.01064vl [cond-mat.mtrl-sci] 6 Jan 2015 


Charge ordering and self-assembled nanostructures in a fee Coulomb lattice gas 
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The compositional ordering of Ag, Pb, Sb, Te ions in (AgSbTe 2 ):c(PbTe) 2 (i_; C ) systems possessing 
a NaCl structure is studied using a Coulomb lattice gas (CLG) model on a face-centered cubic (fee) 
lattice and Monte Carlo simulations. Our results show different possible microstructural orderings. 

Ordered superlattice structures formed out of AgSbTe 2 layers separated by Pb 2 Te 2 layers are ob¬ 
served for a large range of x values. For x = 0.5, we see an array of tubular structures formed by 
AgSbTe 2 and Pb 2 Te 2 blocks. For x = 1, AgSbTe 2 has a body-centered tetragonal (bet) structure 
which is in agreement with previous Monte Carlo simulation results for restricted primitive model 
(RPM) at closed packed density. The phase diagram of this frustrated CLG system is discussed. 

PACS numbers: 64.60.Cn, 81.30.Bx, 81.16.Dn 


I. INTRODUCTION 


Lattice gas with long-range Coulomb interaction 
has attracted considerable interest over the past 10 
years. Two types of long-range models have been 
studied. One where the interaction between the 
charges oc 1/r (Coulomb lattice gas, or CLG), and the 
other where the interaction oc lnr (lattice Coulomb 
gas, or LCG). Studies of various models of one44- 
and two44-dimensional CLG and LCG using differ¬ 
ent methods have shown the existence of multiple 
phase transitions, complexity in phase diagrams and 
their practical applications to real materials, e.g., 
KCu7_ x S 4 r r - Ni 1 ^ x Al x (0H) 2 (C0 3 ) x /2.yH 2 0,... 4 In 

three-dimensional CLG on a simple cubic (sc) lattice, sev¬ 
eral works have been done using either theoretical calcu¬ 
lations (mean-field approximation^ and Pade expansion^) 
or Monte Carlo (MC) simulations^ 9 However, to the 
best of our knowledge, there are no extensive studies of 
CLG on a fee lattice excepting when all the lattice sites 
are occupied by either a positive or a negative charged It 
is well known that fee lattice involves frustration^ Since 
the role played by frustration in the nature of phase tran¬ 
sition in Ising-type systems (on triangular or fee lattice) 
has been of great interest in statistical physics^ - — it is 
of equal interest to see what role frustration effects play 
in long-range Coulomb systems. 


ions is not clear. However, as pointed out by Bile et al .^4 
the electronic structure of these compounds depends sen¬ 
sitively on the nature of structural arrangements of Ag 
and Sb ions. Hence a simple but accurate theoretical 
model is necessary to understand and predict the order¬ 
ing of the ions in these systems. In this paper we present 
a simple ionic model of Ag n Pb m Sb n Te m + 2rl that explic¬ 
itly includes the long-range Coulomb interaction and in 
which the ions are located at the sites of a fee lattice. 
As will be shown in the next section, this problem maps 
onto a spin-1 Ising model on a fee lattice with long-range 
antiferromagnetic interaction. We present details of the 
model in Sec. EH In Sec. m we discuss our Monte Carlo 
simulation results including a full phase diagram in the 
x — T plane. The summary is presented in Sec. IIVI 


II. MODEL 

We use a model where the minimization of electrostatic 
interaction between different ions in the compounds can 
lead to the compositional ordering that exists in the 
system4 7 The total electrostatic energy is then expressed 
as 
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From materials perspective, a quaternary compound 
Ag„Pb m Sb„Te m + 2n has recently emerged as a material 
for potential use in efficient thermoelectric power genera¬ 
tion. It has been found that for low concentrations of Ag, 
Sb and when doped appropriately, this system exhibits a 
high thermoelectric figure of merit ZT of f«2.2 at 800 I<4£ 
It is one of the best known bulk thermoelectrics at high 
temperatures. Quantitative understanding of its proper¬ 
ties requires understanding of atomic structure. Experi¬ 
mental dat a 17 ’ 18 suggest that this system belongs to an 
entire family of compounds, which are compositionally 
complex yet they possess the simple cubic NaCl structure 
on average, but the detailed ordering of Ag, Pb and Sb 


where e is the static dielectric constant, R; T and Qi r 
are, respectively, the position and charge of an atom at 
site t of cell l. This model has been successfully ap¬ 
plied to cubic perovskite alloys4 Here we consider su¬ 
percells of the NaCl-type structure made of two inter¬ 
penetrating fee lattices with possible mixtures of differ¬ 
ent atomic species on Na sites, i.e., r = {Na(Ag,Sb,Pb), 
Cl(Te)}, with periodic boundary condition. Alloying oc¬ 
curs on the Na sublattice. In a simple ionic model of 
Ag n Pb m Sb„Te m + 2n , we can assume the Pb ion to be 
2 + , Te ion to be 2“, Ag ion to be 1 + , and Sb ion to 
be 3 + , i.e., Qi t = {Qz,Naj Qi,c\} = {+1, +3, +2; — 2}, 
where Qi,c\ = 9ci = —2 is independent of l. Focus- 
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ing on the Na sublattice sites where ordering occurs, 
we write Q; TiNa = ^Na + Ag ; , where q Na = +2 and 
A qi = {—1,+1,0}. Substituting the expression for Q;,Na 
into Eq. (1), we can write E = Eq + Ei + E 2 , where 
the subscripts refer to the number of powers of A qi ap¬ 
pearing in that term. Then Eo is just a constant, it is 
the energy of an ideal PbTe lattice; E\ vanishes due to 
charge neutrality. The only term which depends on the 
charge configuration is E 2 ', it is given by 
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where ion positions are measured in unit of the fee lat¬ 
tice constant a ; l and l run over the N sites of the Na- 
sublattice of NaCl structure. Thus if we start from a 
PbTe lattice as a reference system and replace two Pb 
ions by one Ag ion and one Sb ion, we map the sys¬ 
tem unto an effective CLG with effective charges -1, 
+1 (of equal amount) and 0; this implies a constraint, 
Aqi = 0. The model therefore maps onto a spin- 
1 Ising model (s/ = 0, ±1) with long-range antiferro- 
magnetic interaction. The short-range version of this 
model [nearest- (n.n.) and next-nearest-neighbor (n.n.n.) 
interaction] a generalization of this model by adding 
a n.n ferromagnetic interaction^ and a continuum ver¬ 
sion of this model that takes into account the finite size 
of the charged particles (RPM or charged hard sphere 
model)“i2i~— have been investigated. Comparison with 
these works will be made in Sec. mu 

Because of the attraction between +1 and -1 charges, 
the Ag and Sb ions tend to come together and form 
clusters or some sort of ordered structures depending on 
the temperature at which these compounds are synthe¬ 
sized and the annealing scheme. The ordering may be 
quite complex compared to the one on a simple cubic 
lattice because of the frustration associated with spins 
on a fee lattice and antiferromagnetic interaction (in the 
Ising model). In our calculations of Ag ra Pb m Sb n Te m _|_ 2 n> 
an equivalent formula, (AgSbTe 2 )a;(PbTe) 2 (i_ x ), is used; 
where x = 2n/(m + 2n) = l/N^^Aqi] = l/iV]C ; s 2 
(0 < x < 1), is the concentration of Ag and Sb in the Pb 
sublattice. 


III. SIMULATION RESULTS 

To study the thermodynamic properties and mi- 
crostructural ordering of the system, we have done canon¬ 
ical ensemble Monte Carlo simulations following the 
usual Metropolis criterion^ using the energy given by 
Eq. i.e., particles interact via site-exclusive (multiple 
occupancy forbidden) Coulomb interaction. In the Ising 
model problem, this corresponds to a fixed magnetization 
simulation. We used Ewald summation^ to handle this 
long-range interaction employing a very fast lookup table 
scheme using Hoshen-Kopelman algorithm.*^ This model 


is parameter-free in the sense that J = e 2 /ea defines a 
characteristic energy. 

A simulation for a fixed concentration x starts at a 
high temperature with an initial random configuration 
followed by gradual cooling. For each temperature T, we 
use 2 x 10 4 sweeps (MC steps per lattice site) to get ther¬ 
mal equilibration followed by 10 5 sweeps for averaging. 
Particles move either via hopping to empty sites or via 
exchange mechanism. The equilibrium configuration at 
a given temperature T is used as the initial configura¬ 
tion for a study at a nearby temperature. We monitored 
different thermodynamic quantities and look at the mi- 
crostructures. The data presented below were obtained 
with system size L = 8 (i.e., 8 fee cells in one direction, 
2048 lattice sites in total) with periodic boundaries. 



FIG. 1: Energy and heat capacity per particle versus tem¬ 
perature for x = 0.75. Phase transitions occur at T = 0.106 
and 0.21 which are first-order and second-order transitions, 
respectively. There is hysteresis associated with the low T 
transition. 

Figured] shows the energy and heat capacity (obtained 
using energy fluctuation) for x = 0.75 where the two en¬ 
ergy curves correspond to slow cooling and slow heat¬ 
ing. We see evidence of two phase transitions, one at 
T = 0.106 and the other at T = 0.21. The heat capacity 
curve shows peaks at the above two T values. The transi¬ 
tion at higher T is continuous and indicates a lattice gas¬ 
liquidlike phase transition. There is no apparent hystere¬ 
sis associated with this transition. The low T transition, 
on the other hand, appears to be first-order. There is an 
energy discontinuity and there is hysteresis, albeit small, 
associated with this transition. For x < 0.5 we see only 
one transition which is first-order (Fig. [2]). This suggests 
that with decreasing x the system changes from under¬ 
going 2 to 1 phase transition. As x decreases from 0.875, 
the high T continuous and low T first-order phase transi¬ 
tions approach each other and the two transitions merge 
at x « 0.5. As x increases from 0.875, the high T con¬ 
tinuous transition changes to first-order. At x = 1, the 
transition is first-order in agreement with previous sim- 
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FIG. 2: Energy and heat capacity per particle versus tem¬ 
perature for x — 0.25. Phase transition occurs at T = 0.08 
which is first-order. There is hysteresis associated with this 
transition. 

ulation results— We also monitored the structure factor 
S( q) for different q values. For several q values, we find 
that 5(q) changes discontinuously at the first-order tran¬ 
sition and smoothly at a continuous transition. 



FIG. 3: Concentration (x) versus temperature (T) phase dia¬ 
gram constructed from the loci of the heat capacity maxima. 
There are first (solid lines) and second (dotted line) order 
transitions with two possible tricritical points (xt , Ti): A at 
Ri (0.500,0.088) and B at « (0.875,0.290). 

To construct the total phase diagram, we have studied 
energy, heat capacity and structure factor as functions 
of temperature for a series 15 values of the concentra¬ 
tion x. Figure [3] shows the phase diagram constructed 
from the loci of specific heat maxima. The lattice gas- 
lattice solid and lattice liquid-lattice solid transitions are 
first-order. In the limit x = 0, the compound is sim¬ 
ply PbTe, the transition occurs at T = 0 since there are 
no charged particles (effective charges of Pb and Te are 


0). For x = 1, the compound is AgSbTe 2 - Our simula¬ 
tions show a strongly first-order transition at T = 0.38 
and no other transition with decreasing T. This strong 
first-order transition is softened by introducing defects 
into the system (by decreasing x from 1). The hysteresis 
associated with this transition becomes smaller with de¬ 
creasing x from 1 and disappears at x ~ 0.875 showing 
a changeover from a first- to a second-order transition. 
Therefore, we have two possible tricritical points (xt, T t ): 
A at « (0.500,0.088) and B at w (0.875,0.290). More 
accurate results on the tricritical points would require 
further careful large-scale simulations for more number 
of x values, and perhaps much larger systems. 

We would now like to compare our results with those 
of previous simulations carried out for lattice RPM. In 
this model, there is a parameter £ = a /a, where a is 
the hard sphere diameter of the charged particles. For 
£ = 1 which is comparable to our model, Dickman and 
StellZ and Panagiotopoulos and Kumar— have a phase 
diagram for a sc lattice that is similar to ours. They 
found a tricritical point at (x t , T t ) ~ (0.4, 0.14) (Ref. 7) 
and (0.48 ±0.02, 0.15 ± 0.01)— It appears that x t values 
for sc and fee lattices are quite close whereas the T t values 
for the fee lattice is about a factor of 0.6 smaller, perhaps 
due to frustration. As regards the second tricritical point 
(B), it is unique to the fee lattice. Dickman and StellZ 
found a high T continuous phase transition (A-transition) 
in a simple cubic lattice as x increased from 0.4 to 0.82. 
The observation of the high T first-order transition in our 
simulations is similar to the one seen in fully-frustrated 
n.n. and n.n.n. Ising model (x = 1) seen by Phani et 
al^L 



FIG. 4: (Color online) Energy as a function of Monte Carlo 
steps for x = 0.75. The system is quenched from T = 0.3 to 
0.24, then from T = 0.24 to 0.18, and so on. 

In Fig. U we plot the energy of a system with x = 0.75 
as a function of Monte Carlo steps as we quench the sys¬ 
tem from T = 0.3 to 0.002 through several intermediate 
values of T. We use more than 4 x 10 6 moves for each 
T without discarding any step for thermal equilibration. 
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The final configuration at a given T is used as the initial 
configuration for the next T. The fluctuation is large at 
high T and getting smaller with lowering T. From one 
T to another, it takes some time (~ 10 3 steps) for the 
system to equilibrate. As one crosses the transition re¬ 
gion, i.e., from T = 0.24 to 0.18 or from 0.12 to 0.06, 
the result shows the existence of possible local minima in 
energy where the system is in metastable states and then 
goes to a stable state with lower energy. More details on 
quenching studies will be reported in another paper— 



*** ►* *»* ** v** ***** v* *** *** **•* •****■ 



FIG. 5: (Color online) A low temperature configuration for 
x = 0.25 [created using XCrySDen (Ref. 30)]. Dark layers 
are for Ag/Sb, grey layers are for Pb; Te sublattice is not 
shown. This typical configuration showing a ordered super¬ 
lattice structure formed out of AgSbTe 2 layers separated by 
Pb 2 Te 2 layers. 

A typical low temperature structure of 
(AgSbTe 2 ) a: (PbTe) 2 (i_a ; ) is a self-assembled nanos¬ 
tructure with layers of AgSbTe 2 arranged in a particular 
fashion in the PbTe bulk as shown in Fig. [5] for the case 
x = 0.25. Four layers of AgSbTe 2 are separated from 
one another by four layers of Pb 2 Te 2 - This domain is 
again separated by a purely PbTe domain formed by 
eight other layers of Pb 2 Te 2 - Along the z-direction (per¬ 
pendicular to the layers), positive charge and negative 
charge arrange consecutively. This indicates a three- 
dimensional long-range order which is clearly a result 
of the long-range Coulomb interaction. Experimentally, 
high-resolution transmission electron microscopy (TEM) 
images indicate inhomogeneities in the microstructure 
of the materials, showing nano-domains of a Ag-Sb-rich 
phase embedded in a PbTe matri x 17 i 18 which appears to 
be consistent with our results. Also electron diffraction 
measurements show clear experimental evidence of long- 
range ordering of Ag and Sb ions in AgPb m SbTe m+ 2 
(m=18)^ 

In addition to the layered superlattice structures seen 
for several x values, we have also discovered a very inter¬ 
esting structure at x = 0.5, in which an array of tubes 
of AgSbTe 2 and Pb 2 Te 2 are arranged in a checker board 
pattern (Fig. [6]). We find that this structure has the 
same energy as the layered structure consisting of alter- 



FIG. 6: (Color online) A low temperature configuration for 
x — 0.5 [created using XCrySDen (Ref. 30)]. Connected balls 
are for Ag/Sb, unconnected balls are for Pb; Te sub-lattice 
is not shown. Checkerboard pattern formed by AgSbTe 2 and 
Pb 2 Te 2 blocks. 


nate layers of AgSbTe 2 and Pb 2 Te 2 (energy per particle 
E = —1.157278). For x = 1, i.e., AgSbTe 2 , the only or¬ 
dered structure is body-centered tetragonal (bet) struc¬ 
ture with a c-parameter which is double that of the NaCl 
subcell, belonging to space group Mm2. The unit cell of 
this structure has eight ions, with every ion being sur¬ 
rounded by eight ions of opposite charge and four of the 
same charge (Te sublattice is not included here). This 
structure is equivalent to the type-III antiferromagnetic 
structure^ which has been found in n.n. and n.n.n. Ising 
model by Phani et al — It has also been seen in the RPM 
by Bresme et al.- 10 in Monte Carlo simulations and by 
Ciach and Stell— within a field-theoretic approach. 

The comparison we made with a system of size L = 4 
shows no appreciable change in the results for the first- 
order transition except the fact that we did not see any 
hysteresis with L = 4. The energy at a given concen¬ 
tration differs by 0.1% — 0.5% from that obtained for 
L = 8. However, for the continuous transition along the 
line joining A and B in Fig. [3] one expects to see the 
usual finite size effects— Most of the earlier simulations 
have been carried out in systems of similar sizes. For ex¬ 
ample, the system size L = 4 was also used by Bresme et 
al— for a CLG in fee lattice. Bellaiche and Vanderbilt^ 
chose L = 6 for their study of cubic perovskite alloys. For 
a CLG in sc lattice, larger size lattices have been used 
because the number of atoms per unit cell is one in this 
case. For example, Panagiotopoulos and Kumar— and 
Dickman and Stell^ chose L = 12 and 16, respectively. 


IV. SUMMARY 

In summary, our phase diagram has shown the distinct 
feature of having two tricritical points for a CLG in fee 
lattice. We have demonstrated that Monte Carlo simu- 
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lation using an ionic model of (AgSbTe 2 )a;(PbTe) 2 (i-x) 
shows different possible microstructural orderings. We 
have found that layered structures formed out of 
AgSbTe 2 layers separated by Pb 2 Te 2 layers are generic 
low temperature structures. In addition to the layered 
structures, we have also discovered tubular structures for 
x = 0.5. For x = 1, a bet structure has been found, in 
agreement with previous simulation results. Structures 
for other values of x are mixtures of those for x = 0, 
0.5, and 1. These results will be discussed in a separate 


paper isa 
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